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Abstract 

In this paper, wireless systems operating under queueing constraints in the form of limitations on the buffer 
violation probabilities are considered. The throughput under such constraints is captured by the effective capacity 
formulation. It is assumed that finite blocklength codes are employed for transmission. Under this assumption, a recent 
result on the channel coding rate in the finite blocklength regime is incorporated into the analysis and the throughput 
achieved with such codes in the presence of queueing constraints and decoding errors is identified. Performance of 
different transmission strategies (e.g., variable-rate, variable-power, and fixed-rate transmissions) is studied. Interactions 
between the throughput, queueing constraints, coding blocklength, decoding error probabilities, and signal-to-noise 
ratio are investigated and several conclusions with important practical implications are drawn. 

Index Terms: Buffer violation probability, coding rate, decoding error probability, effective capacity, fading 
channels, finite blocklength regime, quality of service constraints, variable-rate/variable-power/fixed-rate transmissions. 
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I. Introduction 

Providing quality of service (QoS) guarantees in the form of limitations on the queueing delays or buffer 
violation probabilities is essential in many delay-sensitive wireless systems, e.g., voice over IP (VoIP), and 
wireless interactive and streaming video applications. Due to the importance of such QoS considerations, 
it is of significant interest to conduct an analysis and provide predictions for the performance levels of 
practical systems. In [1J, effective capacity is proposed as a metric that can be employed to measure the 
performance in the presence of statistical QoS limitations. Effective capacity formulation uses the large 
deviations theory and incorporates the statistical QoS constraints by capturing the rate of decay of the buffer 
occupancy probability for large queue lengths. Hence, effective capacity can be regarded as the maximum 
throughput of a system operating under limitations on the buffer violation probability. 

Recently, there has been much interest in the analysis of the effective capacity of fading channels 
(see e.g., - Ifl2l0 in order to identify the performance of wireless systems operating under statistical 
queueing constraints. However, in almost all prior studies, the service rates of the queueing model (or 
equivalently the instantaneous transmission rates over the wireless channel) are assumed to be equal to 
the instantaneous capacity values although channel coding is performed using a finite block of symbols. 
Moreover, transmissions are assumed to be reliable with no decoding errors. However, it is important to 
note that error-free communication at the rate of channel capacity is generally attained as the codeword 
length increases without bound. Therefore, when finite blocklength codes are employed, transmission is 
necessarily performed in the presence of decoding errors and possibly at rates less than the channel capacity 
in order to have high reliability or equivalently low error probability. 

In [fT3l and Ifl4l . Negi and Goel addressed these considerations. They studied queueing and coding jointly 
and took explicitly into account decoding errors by considering the random coding exponents of error 
probabilities for rates less than the instantaneous channel capacity. For instance, in lfT3l . they analyzed 
the maximization of the joint exponent of the decoding error and delay violation probability through the 
appropriate choice of the transmission rate for given delay bound and constant arrival rate. 

In this paper, we also depart from the idealistic assumptions of communicating arbitrarily reliably at 
channel capacity but follow an approach different from that of [[131 and lfl4l . We consider channel coding 
rates achievable with finite blocklength codes, and incorporate the decoding error probabilities and possible 



1 



retransmission scenarios into the effective capacity formulation. This analysis is facilitated mainly by the 
recent results of Polyanskiy, Poor, and Verdu in [fT51 where the authors identified an approximate maximal 
achievable rate expression for a given error probability in the finite blocklength regime. This expression 
can be regarded as a second-order asymptotic approximation of the channel coding rate at large but finite 
blocklength values. We note that [Tf6l and ifTTl also studied channel coding and achievable error probabilities 
at finite blocklengths by analyzing the mutual information density and its statistics. In iTTvTl . an outage analysis 
is performed by using the distribution of the mutual information density. In |fT8l , a similar outage formulation 
is used to determine the optimal physical-layer reliability and to identify the maximum ARQ throughput. 
On the other hand, neither of the above-mentioned papers have investigated the throughput in the finite 
blocklength regime when the systems operate under buffer constraints. 

Our contributions in this paper can be summarized as follows. We first determine the effective throughput 
in the finite blocklength regime under constraints on the buffer violation probability. Subsequently, we 
study the performance of different transmission strategies. Initially, we consider a scenario in which the 
transmission rate is varied with the fading realizations while the error probability is kept fixed. The optimal 
error probability that maximizes the throughput is shown to be unique. We analyze the impact of the 
power adaptations. Then, we investigate the case in which transmission rate is fixed and error probability 
varies over different transmission blocks. Through numerical results, we analyze the interactions between 
the throughput, queueing constraints, error probabilities, blocklength, signal-to-noise ratio, and different 
transmission strategies. 

The remainder of the paper is organized as follows. Section [II] describes the fading channel model. In 
Section Unl we provide preliminaries on the effective capacity as a measure of the throughput under statistical 
QoS constraints. In Section [IVl we provide our results on the effective throughput in the finite blocklength 
regime. We conclude in Section |Vj Several proofs are relegated to the Appendix. 

II. Channel Model 

We consider a frequency-flat channel model, and assume that the fading coefficients stay fixed for a block 
of m symbols and then change independently for the following block. Under this block-fading assumption, 
the channel input-output relation in one coherence block can be expressed as 
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y = hx + n 



(1) 



where x and y are the m-dimensional, complex, channel input and output vectors, respectively. The input 
is subject to an average power constraint, i.e., E{||x|| 2 } < mP . h is the complex- valued fading coefficient 
with finite second moment, i.e., E{|/i| 2 } < oo. We assume that both the receiver and transmitter have 
perfect channel side information (CSI) and hence perfectly know the instantaneous realizations of the fading 
coefficients. However, the assumption of perfect CSI at the transmitter is relaxed in Section ITV-B1 Finally, n 
represents the Gaussian noise vector whose components are independent and identically distributed (i.i.d.), 
complex, circularly symmetric, Gaussian random variables with mean zero and variance N , i.e., n ~ 
£/V(0, iV I m ) where I m denotes the m x m identity matrix. 

III. Throughput under Statistical Queueing Constraints 

In HJ, Wu and Negi defined the effective capacity as the maximum constant arrival rate that a given service 
process can support in order to guarantee a statistical QoS requirement specified by the QoS exponent 9 P. 
If we define Q as the stationary queue length, then 9 is the decay rate of the tail of the distribution of the 
queue length Q: 

Urn = (2) 

q^oo q 

Therefore, for large q max , we have the following approximation for the buffer violation probability: 

Hence, while larger 9 corresponds to more strict QoS constraints, smaller 9 implies looser QoS guarantees. 
Similarly, if D denotes the steady-state delay experienced in the buffer, then P(D > rf max ) ~ e~ e<Wmax for 
large d max , where S is determined by the arrival and service processes 0. Therefore, effective capacity 
formulation provides the maximum constant arrival rates that can be supported by the time- varying wireless 
channel under the queue length constraint P(Q > g max ) < e - eqmax for large q max or the delay constraint 
P(D > (i m ax) < e ~ 6,<Wmax for large d max . Since the average arrival rate is equal to the average departure rate 

'For time- varying arrival rates, effective capacity specifies the effective bandwidth of the arrival process that can be supported by the channel. 
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when the queue is in steady-state [4], effective capacity can also be seen as the maximum throughput in the 
presence of such constraints. 

The effective capacity is given by (Q]|, CD, 0) 

R E = - \im Uog e E{e- es ^} (3) 

t-s>oo tit 

where S[t] = J2i=i^i * s me time-accumulated service process and {Ri,i = 1,2, . . .} denotes the discrete- 
time stationary and ergodic stochastic service process. We would like to note that in the remainder of the 
paper, we will refer to Re as the effective rate rather than the effective capacity since Re in our setup is the 
throughput when the service rates are equal to the approximate channel coding rates in the finite blocklength 
regime. 

IV. Effective Throughput with Finite Blocklength Codes 

In |fT5l . the authors have studied the channel coding rate in the finite blocklength regime. For general 
classes of channels, they have obtained new achievability and converse bounds on the coding rate for a given 
finite blocklength and error probability. In particular, for the real, additive white Gaussian noise (AWGN) 
channel, the transmission rate (in bits per m channel uses) with error probability < e < 1, signal-to-noise 
ratio (snr), and coding blocklength m is shown to have the following asymptotic expression lfl"5l Theorem 
54]: 

7TI / 777 / I \ 

r = - log 2 (l + snr) - J- ( 1 - (sNR+1)2 j Q-\e) log 2 e + O(logm) (4) 

where Q(x) = -j= e~* 2//2 dt is the Gaussian Q-function. Denoting the rate in bits per channel use by f, 
we can write 



L = Il 0g2 (l + snr) - J± (l - -4^V 'M log 2 e + (5) 
m 2 v 2m \ (sNR+lr/ 



m 



« \ log 2 (l + snr) -J±(l- 1 ) Q~\e) log 2 e (6) 
2 y 2m \ (snr + l) z J 

where the approximation is accurate for sufficiently large m. Note that the above results are for the AWGN 
channel with real input and real output. 

In this paper, we consider a fading Gaussian channel model with complex- valued input and output, and 
assume that channel coding is performed in each coherence interval of m symbols, during which the fading 
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stays fixed. Under these assumptions, coding over a fading Gaussian channel can be seen as coding over 
a real Gaussian channel (with a certain channel gain) using a coding blocklength of 2m. The following 
arguments provide a detailed description of this approach. Knowing the channel fading coefficient h, the 
receiver can multiply the received signal with e~^ h , where 9h is the phase of h, and obtaiij^l 

y = yr+jy< = ye~ jdh = |ft|x + n = |/i|x r + h r + j(\h\xi + fi<) (7) 

where y r , x r , n r and fi, Xj, fij denote the real and imaginary components, respectively, of the output vector 
y, input vector x, and noise vector n. It can be easily verified that n = ne~^ 9h has the same statistics as n 
and hence n ~ CJ\f(Q, N l m ). Now, the above channel input-output relation can also be written as 

fi] = \H [ X r Xj] + [Sp fii] (8) 

where [y r fi] denotes the vector formed by concatenating yy and yV Since the real and imaginary components 
are m-dimensional vectors, the above channel model is a real Gaussian channel with 2m dimensional input 
and output and with channel gain \h\. Note that the real and imaginary noise components n r and are 
independent due to the assumption of the circular symmetry of the additive complex Gaussian noise. For this 
channel, the coding rate (in bits per m channel uses) in the i th block achieved with block error probability 
e is 

n =mlog 2 (l + SNR|^| 2 ) - (l - 7j ^| fe ,| 2 + 1 )2 ) < ^" 1 ( £ ) log 2 e + 0(\og2m) (9) 

where hi denotes the fading coefficient in the i th block. Note that the expression in © is obtained from that 
in © by replacing m with 2m, and snr with SNR|/ij| 2 = -^-|/ij| 2 , which is the received signal-to-noise ratio 
in the i th block. Now, the normalized rate in bits per channel use is approximately 



rn 



f i = ^ = log 2 (l + snr|^| 2 ) - J- (l- ) W(e)log 2 e (10) 

U m \ (SNR|/lj| z + \y J 



for large enough m for which °^ lo ^ 2m ' ) is negligible. Henceforth, we assume that the instantaneous trans- 
mission rate in each coherence block of the fading channel is given by the expression in (flOl) . Since the 



2 Note that multiplication of the channel output with the e jBh just rotates the output, is a reversible operation, and hence does not lead to 
any loss of information. 
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block error rate is e, this rate is attained with probability 1 — e. We assume that the receiver reliably 
detects the errors, employs a simple ARQ mechanism and sends a negative acknowledgement requesting the 
retransmission of the message in case of an erroneous reception. Therefore, the data rate is effectively zero 
when error occurs. Under this assumption, the service rate (in bits per m channel uses) in each block is 

{0 with prob. e 
(11) 
mvi with prob. (1 — e) 

With the above service rate characterization, we immediately obtain the following expression for the effective 
rate. 

Proposition 1: The effective rate (in bits per channel use) at a given snr, error probability e, blocklength 
m, and QoS exponent 9 is 

Re{0) = -^jlog e E w2 {e+ (1 - e)e- e ™} (12) 

where f is given in (flOl) and the expectation is with respect to \h\ 2 . 

Proof: We first note that the service rate {Ri} is an i.i.d. process due to the facts that the fading process 
is i.i.d. in different blocks and the noise is an i.i.d. process leading to the independence of error events in 
different blocks. Now, we have 

R E {6) = - lim ilog e E{e- es M} (13) 
i->oo yt 

= - lim — log e E{e- e Rl } (14) 

=-,'^s ioE - E {n^ H '} < i5 > 

i * 

= -lim-log e TTE{e-^} (16) 

i=l 

= -limitog e (E{e-^}) t (17) 

= -lim^-tlog e E{e-^} (18) 

t-s>oo yt 

= -ilog e E{e-^} (19) 

= -ilog e E W2 { e +(l- e )e^} (20) 
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Above, (fT6l) follows from the independence of the service process and (TT71 ) is due to its being identically 
distributed. The expression inside the expectation in (1201 is obtained by evaluating the expected value of 
e -0Ri f or fj xe( j |/^|2 Finally, (fT2)) is obtained by normalizing (|20~1) by m to have the effective rate in the units 
of bits per channel use, and by dropping the time index i. □ 

Note that the effective rate is a function of the QoS exponent 8, blocklength m, signal-to-noise ratio 
snr and error probability e. Since we assume that coding is performed in each coherence interval, the 
blocklength m is determined by the statistics of the fading process. The value of 8 can be dictated by the 
application requirements and snr depends on the power budget. Given the values of these parameters, the 
remaining parameter e can be optimized to maximize the throughput. Note that large e implies that the 
transmitter attempts to transmit the data at a high rate but at the risk of more frequent errors and hence 
retransmissions. On the other hand, if e is small, the instantaneous transmission rate is low but the reliability 
of the transmissions is high. The following result shows that the optimal e is unique. 

Proposition 2: Assume that the values of m, 8 > 0, and snr > are fixed. Then, the function 

M/(e) = E|,| 2 {e + (l-e)e- emf } (21) 

is strictly convex in e and therefore the optimal value of e that minimizes this function or equivalently 
maximizes the effective rate in (fT2l) is unique. 
Proof: See Appendix [A] 

Note that the convexity result indicates that the optimal error probability e* can be easily found using 
standard convex optimization methods. The analysis and the resulting e* provide guidelines on the design 
of the channel codes and their strength. Note further that the above result is shown for the case in which 
8 > 0. If there are no QoS constraints and hence = 0, then we have the following corollary to Proposition 
I 

Corollary 1: When 8 = 0, the effective capacity becomes 

R E (0) = lim R E {8) = (1 - e)E N2 {f } (22) 

where f is given in (flOi 

Note that the Re(0) is the average transmission rate averaged over the fading states. Below, we show that 
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Re(0) is a strictly concave function of e. 

Proposition 3: Assume that the values of m, and snr > are fixed. Then, the function 

R E (0) = (1 - e)E|„| 3 {f} (23) 

is strictly concave in e and therefore the optimal value of e that maximizes this effective rate is unique. 
Proof: See Appendix [Bj 

Next, we provide numerical examples to illustrate the results. Although the preceding analysis is applicable 
to any fading distribution with finite power, we consider a Rayleigh fading channel in the numerical analysis, 
and assume that the fading power z = \h\ 2 is exponentially distributed with unit mean (i.e., has the probability 
density function f z (z) = e~ z ). 

In Figure \T\ we plot ^(e) = E^a {e+ (1 — e)e _6,mr } as a function of the error probability e in the 
Rayleigh fading channel. In the figure, snr = dB and the blocklength m = 1000. We provide curves for 
different values of the QoS exponent 9 > 0. In all cases, we immediately observe the strict convexity of the 
curves, confirming the result in Proposition |2] Indeed, the optimal error probabilities that minimize \&(e) are 
unique and are equal to e* = 0.0127, 0.0061, 0.0084 for 9 = 0.001, 0.01, 0.1, respectively. 

In Fig. [2l we plot the effective rate in (fl~2l) as a function of the error probability e. The other parameters 
are the same as in Fig. \T\ Notice that we have also included in this figure the throughput curve for the case 
in which 9 = 0. Note that if 9 = 0, the system does not have any queueing constraints. In Proposition |3j we 
have shown that Re(0) is a strictly concave function of e and the optimal e* that maximizes Re(0) is unique. 
The strict concavity is observed in Fig. [2] The optimal value of the error probability in the case of 9 = is 
e* = 0.0171. For 9 > 0, the effective rate curves are not necessarily concave. In Fig. |2l we observe that these 
curves are quasiconcave and, as predicted by Proposition [H they are maximized at a unique e*. The optimal 
error probabilities for the cases in which 9 > are equal to the same ones obtained in Fig. \T\ At the optimal 
error probabilities, the maximum effective rate values are Re = 0.7750, 0.6256, 0.2246, 0.0329 bits/channel 
use for 9 = 0,0.001,0.01,0.1, respectively. Note that increasing 9 leads to more stringent QoS constraints, 
and we observe that the effective rate and hence the effective throughput diminishes as 9 increases. This 
trend is also clearly seen in Fig. [3] where we plot the maximum effective rate values (i.e., effective rate at 
the optimal error probability e*) as a function of 9. 
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Another interesting analysis is the behavior of e* as a function of 9. This is depicted in Fig. HI Here, 
we observe that as 9 increases and therefore the QoS limitations become more stringent, the value of e* 
initially decreases sharply. Hence, the transmitter opts for more reliable but low-rate transmissions. On the 
other hand, as 9 increases beyond approximately 0.028, the trend reverses and e* starts to increase. The 
transmitter increases the transmission rate at the cost of increased e* and hence more retransmissions. When 
9 exceeds 0.298, e* starts decreasing again. Note that for high values of 9, the effective rate is small. This 
small effective rate can be supported by low-rate transmissions. Hence, when 9 is high beyond a threshold, the 
transmitter chooses to transmit at low rates and keep the error probability and the number of retransmissions 
low as well. 

In Fig. |5l we plot the effective rate as a function of the blocklength m for 9 = and 9 = 0.001. The 
solid-lined curves correspond to the effective rate in (PT21) optimized over e. The dashed curves correspond 
to the effective rate of the ideal model in which the service rate is equal to the instantaneous capacity, i.e., 

f = log 2 (l + snr|/i| 2 ), (24) 

and the error probability is assumed to be zero, i.e., e = 0. Here, we have interesting observations. When 
9 = and the ideal model is considered, then the effective rate is Re(0) = E| /l |2{log 2 (l + snr|/i| 2 )}, which 
is the ergodic capacity of the fading channel and is clearly independent of the blocklength. On the other 
hand, if the service rate is given by f in (flOl ), the effective rate -Re(O) = (1 — e)E| /l |2 {r} increases with 
blocklength m as seen in Fig. |5J In the presence of QoS constraints, i.e., when 9 > 0, we have stark 
differences. Under the idealistic assumption of transmitting at the instantaneous capacity with no errors, we 
see from the behavior of the dashed curve for 9 = 0.001 that effective rate decreases with increasing m. 
The reason is that since m is the coherence duration over which the fading state remains fixed, larger m 
corresponds to slower fading and slow fading is detrimental for buffer-constrained systems. In a slow-fading 
scenario, deep-fading can be persistent causing long durations of low rate transmissions leading to buffer 
overflows. In the finite blocklength regime, as seen in the behavior of the solid-lined curve of the case of 
9 = 0.001, there is a certain tradeoff. Initially, increasing m improves the performance as this allows the 
system to perform transmissions with longer codewords and to have higher transmission rates. However, if 
m increases beyond a threshold, slowness of the fading starts to degrade the performance. 
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In all cases in Fig. [5J the gap between the dashed and solid-lined curves diminishes as m increases 
since the idealistic model becomes more accurate. On the other hand, for moderate values of m (e.g., when 
m < 2000), the idealistic assumptions lead to significant overestimations of the performance. 

Finally, we provide numerical results for the optimal effective rate and optimal error probability as a 
function of snr in Figs. [6] and [7J respectively, for 9 = 0,0.001, and 0.01. We see that, for fixed 9, 
increasing the snr improves the throughput and also the reliability of the transmissions by lowering the 
error probabilities. 

A. The Impact of Power Adaptation 

Heretofore, we have considered the scenario where the transmitter knows the fading coefficients {hi} and 
performs variable-rate transmission with the same average power P in each coherence block of m channel 
uses. In this section, we investigate the gains achieved by varying the transmission power as well with 
respect to fading. Let us denote the power adaptation normalized by the noise power by //(snr, 9, \h\ 2 ). 
With this adaptation policy, the transmission rate is 

f = load + //(snr, B, WW - ^ ZZZZ^ Z^Q-^ log2e (25) 

which is obtained by replacing snr with //(snr, 9, \h\ 2 ) in (TTOb . Finding the optimal power adaptation policy 
that maximizes f or the effective rate Re{9) = —-^ log e E^p {e + (1 — e)e _(9mr } is in general a difficult 
task due to the facts that both the first and second terms on the right-hand side of (l25l) are concave functions. 
Hence, f is neither concave or convex. For this reason, we resort to suboptimal strategies. One viable policy, 
//*, is the one that maximizes the effective rate when the service process is assumed to be equal to the 
instantaneous capacity log(l + //(snr, 9, \h\ 2 )\h\ 2 ) with zero error probability, i.e., 

//*(snr,0, \h\ 2 )} = arg max -±l g e E, h , a J e -^i° g2 (i+MSNRAN 2 )N 2 )) . (26) 

//* is derived in [5] and is given by 

f— — 1 — l^l 2 >« 

a j+i { \ h \2)j+i ' ' (27 ) 
\h\ 2 <a 
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where (5 



dm 
log„2 



and a is chosen such that the average long-term signal-to-noise ratio constraint, 



S| fe i2{/i(sNR, 9, \h\ 2 )} < snr is satisfied with equality. Note that this policy is close to the optimal one when 
the blocklength is large and hence f is close to log(l + //(snr, 9, \h\ 2 )\h\ 2 ) and e is close to zero. 

In Fig. [H the optimal effective rate is plotted as a function of 9 for both fixed- and variable-power cases. 
In the fixed power case, snr = dB in each coherence block. When power adaptation is employed, signal- 
to-noise ratio /a (snr, 9, \h\ 2 ) varies in each block while satisfying E^pj/^SNR, 9, \h\ 2 )} < snr = dB. The 
improved performance with power control is observed in the figure. 



B. Fixed-Rate Transmissions 

The analysis above has assumed that the transmitter has perfect knowledge of the fading coefficients and 
can perform variable-rate and/or variable-power transmissions in each coherence block. On the other hand, 
it is practically interesting to consider cases in which the transmitter does not know the channel and send the 
information at a fixed rate. Additionally, the transmitter may prefer fixed-rate transmissions, even when it 
knows the channel, due to complexities in varying the transmission rate for each block. Motivated by these 
considerations, we assume in this section that the transmitter sends the information at the fixed rate f/. 
Under this assumption, error probability e varies with the fading realizations. The analysis in the previous 
sections have, on the other hand, considered the scenarios in which the error probability is fixed for all 
channel states. 

From (flOl ), which provides the fundamental tradeoff between the rate and error probability in the finite 
blocklength regime, we can easily see that the error probability for fixed f/ is 



Q 



( \ 

log 2 (l + SNr|/?,| 2 ) — Tf 



1 



(SNR|h| 2 +l) 2 



log 2 e 



(28) 



Note that e is a function of the fading magnitude \h\, signal-to-noise ratio snr, and blocklength m. The 
service rate (in bits per m channel uses) is now 



Ri = < 



with prob. e = Q 
mff with prob. (1 — e) 




log 2 (l+SNR\h\ 2 )-r f 



(SNR|h 



^w) log2e 



(29) 



11 



It can also be immediately seen that for given snr, blocklength m, QoS exponent 9, and fixed-rate r/, the 
the effective rate in bits per channel use is 



Re{0) 



1710 



log e E N2 { e + (l-6)e- emf /} 



(30) 



which is essentially the same as in (|12)) . The only difference is that we now have the rate fixed and error 
probability varying. Similarly, when = 0, we have 



/ 



R E (0) = E| fc|2 {(l - e)f f } = (1 - E W 2{e})f f 



1 - E|M2 I Q 



( \\\ 

log 2 (l + SNR|/l| 2 ) — Tf 



\ 



m ( (SNR /, -+] i- 



It is instructive to investigate what is obtained as m — > oo. We immediately see that 



lim Q 



( \ 

log 2 (l + SNR|/i| 2 ) — Tf 



leading 



A 



i - 



(SNR|/i| 2 +l) 2 



log 2 e 



if r f < log 2 (l + snr|/i| 2 ) 

1 if ff > log 2 (l + snr|/i| 2 ) 



lim E| h | 2 { Q 



( \ 

log 2 (l + SNR|/l| 2 ) — Tf 



(SNR|/i|2+T)^ 



log 2 e 



P(f f >log 2 (l + SNR|/ i | 2 )) =P out . 



Therefore, in the limit as m — > oo, 



/ 



r f- 



(31) 



(32) 



(33) 



R E (0) ->• (1 - Pout)f/ 



(34) 



which is defined as the capacity with outage 11201 Section 4.2.3]. Therefore, Re(0) in (|3TI) can be seen as 
the outage capacity in the finite blocklength regime. Furthermore, Re(0) in (|3Q|) can be regarded as the 
generalization of such a throughput measure to the scenario with QoS limitations. 

In Figs. H] - [HI we illustrate the numerical results. In Fig. effective rate is given as a function of 



The interchange of the limit and the integral (or equivalently the expectation) can be easily justified by noting the boundedness of the 
Q-function, i.e., \Q(-)\ < 1, and invoking the Dominated Convergence Theorem. Additionally, we implicitly assume that the random variable 
log 2 (l + SNR|/i| 2 ) does not have a mass at ff and hence ff — log 2 (l + SNR|h,| 2 ) is a zero-probability event and this event does not affect the 
expectation. 
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the fixed transmission rate ff. We observe that the effective rate curves are quasiconcave and moreover 
they are maximized at a unique value of ff. We also observe that the maximum value of the effective rate 
diminishes with increasing 9. This is more clearly seen in Fig. [TOl where the optimal effective rates (optimized 
over ff) are plotted as a function of 9. In this figure, we have curves for both fixed-rate and variable-rate 
transmissions. Effective rate for the variable-rate transmission is computed by maximizing (fT2l) over e. It is 
interesting to observe that fixed-rate transmissions perform worse than variable-rate transmissions for small 
values of 9. However, for 9 > 0.13, fixed-rate transmissions start outperforming. Hence, for high enough 
values of 9, fixing the transmission rate and having the error probability vary in each block provide better 
performance than requiring the error probability to be fixed by varying the rate. Finally, in Fig. [TT] we note 
that as 9 increases, the optimal fixed rate ff, which maximizes Re(9) in (1301 , diminishes. 

C. Sending Independent Messages over Two Parallel Channels 

So far, we have assumed that the transmitter sends a single codeword x = [x r Xj] of length 2m in 
m channel uses. Another approach is to transmit two independent messages using codewords x r and Xj 
selected from two independent codebooks. Note that now the codeword length is m. These two independent 
codewords can be seen to be sent through two independent parallel channels: 

y = fr+jfi = ye~ 36h = |/i|x + n = J h | x r . + h r +j ( | h | Xj + n» ) ■ (35) 

channel 1 channel 2 

Since the blocklength is m for each codeword, the transmitter sends the information through each channel 
in the i th block duration at the following rate with block error probability e: 



n, P =\ log 2 (l + snrN 2 ) - (l - (SNR| J 2 + 1)2 )^ 1 (^) log 2 e (36) 

where the subscript p is introduced to differentiate this rate from that in (TTOT ). Since errors occur independently 
in each channel, the service rate (in bits per m channel uses) in each block duration of m channel uses is 



Ri 



with prob. e 2 

mfjp with prob. 2e(l — e) • (37) 
2mf i)P with prob. (1 — e) 2 
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Effective rate for this service rate can easily be found as in the proof of Proposition [T] and the proof of the 
following result is omitted for brevity. 

Proposition 4: When the transmitter sends two independent messages over the independent real and 
imaginary channels, the effective rate in bits per channel use at a given snr, error probability e, blocklength 
m, and QoS exponent 9 is 

Re(0) = ~ hg e E ]h]2 {e 2 + 2e(l - e)e~ 6m ^ + (1 - e) 2 e" 2 ^)} (38) 
= -^0 lo ^ E ^{( e+ ^-^ emfp ) 2 } (39) 

where r p is given in (|36l) . 

In this case, it can again be easily shown that the error probability e that maximizes the effective rate in 
(|39| ) is unique. The following is a corollary to Proposition [2l 

Corollary 2: Assume that the values of m, 9 > 0, and snr > are fixed. Then, the function 

%(e) = E w > {(e + (1 - e)e- em ^) 2 } (40) 

is strictly convex in e and therefore the optimal value of e that minimizes this function or equivalently 
maximizes the effective rate in (|39l ) is unique. 
Proof: See Appendix O 

In the absence of QoS constraints, the effective rate becomes 

R E (0) = (1 - e)E W2 {2f p } (41) 



- (1 - e)IW |log 2 (l + SNR |^| 2 ) - sjl (l " (^^ ) Q " 1(£)1 ° g2e ) ' (42) 

which can immediately be seen to be smaller than the effective rate in (f22|) . Hence, when 9 = 0, using 
two codewords, each of length m, provides lower throughput than using a single codeword of length 2m. 
Surprisingly, as we observe in Fig. [121 the throughput achieved by sending two codewords is higher if 
9 increases beyond a threshold. Therefore, under strict QoS constraints, sending in each coherence block 
multiple codewords with shorter lengths may be preferable. 
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V. Conclusion 



We have analyzed the performance of buffer-constrained wireless systems in the practical scenario in which 
transmissions are performed using finite blocklength codes with possible decoding errors at the receiver. 
Employing a recent result on coding rate in the finite blocklength regime, we have determined the effective 
rate expression as a function of the QoS exponent, coding blocklength, decoding error probability, and 
signal-to-noise ratio, and characterized the throughput under statistical QoS constraints. We have discussed 
different transmission strategies. In the case in which the transmission rate is varied and the error probability 
is kept fixed across different fading realizations, we have shown that the effective rate is maximized at a 
unique error probability. This optimal decoding error probability gives us insight on the required reliability 
of the channel codes. Through numerical results, we have investigated how the optimal effective rate and 
optimal error probability vary with the QoS exponent 9. We have also had interesting observations on the 
performance as a function of the blocklength. We have analyzed the throughput improvements through 
power adaptation. We have studied the practical scenario in which the transmitter sends the information at 
a fixed-transmission rate. We have seen that while variable-rate schemes provide higher effective rate at 
low values of 9, fixed-rate transmissions start performing better as 9 increases. Finally, we have noted that 
sending multiple codewords with shorter blocklengths in each coherence interval can become a favorable 
strategy under stringent QoS constraints. 



Appendix 



A. Proof of Proposition [2] 



We first prove the following Lemma. 



Lemma 1: For fixed m, snr > 0, and \h\ 2 > 0, 



f(e) = (1 - e)e 



— 8mr 



(43) 



is a strictly convex function of e. 



Proof: We first express 



9mr = aQ~ l (e) + b 



(44) 
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where, from (flOl) , 



a = 9\ m ( 1 — 



1 



(snrI/iP + I) 2 



loge and b = —9mlog 2 (l + sm\h\* 



(45) 



Note that since snr > 0, \h\ 2 > and 9 > 0, we have a > 0. With the above definitions, we can write 

/( e ) = (1 - e )e aQ ' 1{t)+h . (46) 
The first and second derivatives of /(e) with respect to e can easily be found as follows: 



/(e) = o(l - e)Q-\e) - 1 



D aQ- 1 (t)+b 



/(e) = a(l - eXQ" 1 ^ - 2Q- 1 (e) + (1 - e )g- l (e) 



ae 



aQ _1 (e)+6 



(47) 
(48) 



where Q _1 (e) and Q^ 1 (e) denote the first and second derivatives, respectively, of Q _1 (e) with respect to 
e. Next, we employ several techniques used in lfT8l Appendix A] to prove the Lemma. Note that for an 
invertible and differentiable function g, we have g(g~ 1 (x)) = x. Taking derivative of both sides of this 
equality leads us to 



g-\x) 



(49) 



g(g 

where g~ l {x) denotes the derivative of g~ l with respect to x, and g(g^ 1 (x)) is the derivative of g evaluated 
at g~ l (x). Following this approach and noting that 



Q(x) 



POO 

/ —=e' t2/2 dt, and Q(x) 

J x V27T 



-x 2 /2 



(50) 



we can easily find the following expression: 



Q _1 (e) = -V2ne 



(51) 



Note that Q x (e) < for any < e < 1. Differentiating Q x (e) with respect to e, we obtain the second 
derivative as follows: 



Q-\e) = 27iQ-\e)e 



(52) 
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Next, we consider two cases: 

1) e < 1/2: First, we assume that e < 1/2. Under this assumption, we have Q~ l (e) > and hence 
Q^ 1 (e) > 0. Together with the fact that Q~ l (e) < 0, we immediately see that 

/(e) > for e < 1/2. (53) 

2) e > 1/2: Next, we analyze the case in which e > 1/2 and therefore Q _1 (e) < 0. We concentrate on 
the term inside the square parentheses in (|48T) . Using (|5TT) and (l52l) . and defining x = Q _1 (e) or equivalently 
Q(x) = e, we can write 

a(l - e){Q-\e)f - 2Q~\e) + (1 - e)Q-\e) (54) 

= a(l - e^vre^ 1 ^ 2 + l^e^^l 2 + (1 - e^Q- 1 ^^ 1 ^)) 2 (55) 

= a(l - g(x))27re x2 + 2v / 27e l2/2 + (1 - Q(x))2vr:re x2 (56) 

= e x2/2 (2n(l - Q{x)){x + a)e x2/2 + 2v^tt) (57) 

> e x2/2 (2ti(1 - Q(x))xe x2/2 + 2y/2^ (58) 

> e - 2 /2 ( 2?r * e - x2 / 2 xe x2 ' 2 + 2y/2^] (59) 



> e x 12 (-v 7 ^ + 2V2^j (60) 

> e x2/2 (v 7 ^) > 0. (61) 



Above, (1581) follows from the fact that a > and hence x + a > x. (1591) is obtained by using the upper 
bound, 

1 - Q(x) = Q(-x) < -=J e~ x2/2 for x < 0, (62) 

V27r(— x) 

and recognizing that by our assumption x = Q~ l (e) < 0, and (1 — Q(x)) is multiplied above by x < 0, 
enabling us to find a lower bound. From the above discussion, we conclude that 

/(e) > for e > 1/2. (63) 
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Finally, note that when e = 1/2 and hence Q 1 (e) = Q 1 (1/2) = 0, we have 

o(l - eXQ- 1 ^)) 2 " 2g- X (e) + (1 - e)^) (64) 
= a(l - e)2ne ( - Q ~ 1 ^ 2 + 2^e^~ 1 ^' 2 + (1 - e)2nQ-\e)e^~ 1 ^ a (65) 
= avr + 2^ > 0, (66) 

and therefore /(1/2) > 0. Since /(e) > for all e G [0, 1], /(e) is a strictly convex function of e. □ 
We now define 

^(e) = e + /(e) = e + (l-e)e- 9mf (67) 

which is also strictly convex as it can be immediately seen that ip(e) = /(e) > for snr > and \h\ 2 > 0. 
Note that if either snr = or \h\ 2 = 0, the coding rate becomes f = 0, leading to ip(e) = 1. Since the 
nonnegative weighted sum of strictly convex functions is strictly convex [191 and since the addition of a 
constant (in the case of \h\ 2 = 0) does not have an impact on the strict convexity, we immediately conclude 
that 

*(e) = E, A | a Me)} = E| h , a {e + (1 - e)e~ dmf } (68) 
is strictly convex in e, proving Proposition [21 ■ 

B. Proof of Proposition \3\ 

The proof is similar to that of Proposition |2] in Appendix [A] and will be kept brief. Let's first consider 
the function 

0(e) = (1 - e)f = (1 - e) ( Cl - c 2 Q-\e)) (69) 



where we define c\ = log 2 (l + snr|/i| 2 ) and c 2 = y ^ ^1 — (snrJ^P+T) 7 ) l°g e - Note that if either snr = 
or \h\ 2 = 0, then c\ = c 2 = and 0(e) = for all e. Next, we consider the case in which snr > and 



18 



\h\ 2 > 0, and therefore c\ > and c 2 > Q The second derivative of 0(e) with respect to e is 

0(e) = 2c 2 Q-\e) + c 2 (e - l)g -1 (e). (70) 

Using similar arguments as in Appendix lAl we can easily see that for e < 1/2, 0(e) < 0. For e > 1/2, we 
can show, employing steps similar to those in (|54|)-(|6T1). that 

0(e) < -c 2 V2^e x2/2 < (71) 

where x = Q~ 1 (e). When e = 1/2, we have 0(e) = —2^/2nc 2 < 0. Since 0(e) < for all e, 0(e) is a 
strictly concave function of e when \h\ 2 > and snr > 0. As argued similarly in Appendix lAl since the 
nonnegative weighted sum of strictly concave functions is strictly concave [fT9l and since the addition of a 
constant (in the case of \h\ 2 = 0) does not have an impact on the strict concavity, we conclude that 

R E (0) = (1 - e)E, h ,,{f} = (1 - e)E, h ,.{(c a - c 2 Q-\e))} = E„| 2 {0(e)} (72) 

is a strictly concave function of e. ■ 

C. Proof of Corollary |2] 

From the proof of Proposition [2] in Appendix lAl it immediately follows that e + (1 — e)e~ dmTp is a strictly 
convex function of e. Then, (e + (1 — e)e~ 9mVp ) 2 is strictly convex due to the facts that f(x) = x 2 is a strictly 
convex and increasing function of x and the composition f(g(x)) is strictly convex function when g(x) is 
a strictly convex function |fT9l Section 3.2.4]. Then, strict convexity of E^p |(e + (1 — e)e~ dmTp ) 2 } follows 
from the arguments employed at the end of Appendix [A] ■ 
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Fig. 1. The function ^(e) vs. the error probability e in the Rayleigh fading channel. SNR = dB and the blocklength is m = 1000. 
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Fig. 2. Effective rate Re vs. the error probability e in the Rayleigh fading channel. SNR = dB and the blocklength is m = 1000. 



21 




Fig. 3. The optimal effective rate Re vs. QoS exponent 9 in the Rayleigh fading channel. SNR = dB and the blocklength is m — 1000. 
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Fig. 4. The optimal error probability e* vs. QoS exponent 8 in the Rayleigh fading channel. SNR = dB and the blocklength is m = 1000. 
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Fig. 5. The optimal effective rate Re vs. the blocklength m in the Rayleigh fading channel. SNR = dB and the QoS exponent is 9 — 0.001. 
Dashed curves correspond to the effective rate of the ideal model in which the service rate is equal to the instantaneous channel capacity and 
error probability is zero. 
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Fig. 9. Effective rate Re vs. the fixed transmission rate 77 in the Rayleigh fading channel. SNR = dB and the blocklength is m — 1000. 




Fig. 10. The optimal effective rate Re vs. 6 in the Rayleigh fading channel for both variable-rate and fixed-rate transmissions. SNR = dB 
and the blocklength is rn — 1000. 
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Fig. 11. The optimal fixed-transmission rate 77 vs. QoS exponent 6 in the Rayleigh fading channel. SNR = dB and the blocklength is 
m = 1000. 




Fig. 12. The optimal effective rate Re vs. 8 in the Rayleigh fading channel. SNR = dB and the blocklength is m = 1000. The dashed 
curve is the effective rate in j 1 2t maximized over e and the solid curve is the effective rate in ( 1391 ) maximized over e. 
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